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We study the effects of Rashba spin-orbit coupling on the spin torque induced by 
spin waves, which are the plane wave dynamics of magnetization. The spin torque 
is derived from linear response theory, and we calculate the dynamic spin torque by 
considering the impurity-ladder-sum vertex corrections. This dynamic spin torque is 
divided into three terms: a damping term, a distortion term, and a correction term 
for the equation of motion. The distorting torque describes a phenomenon unique 
to the Rashba spin-orbit coupling system, where the distorted motion of magnetiza¬ 
tion precession is subjected to the anisotropic force from the Rashba coupling. The 
oscillation mode of the precession exhibits an elliptical trajectory, and the ellipticity 
depends on the strength of the nesting effects, which could be reduced by decreasing 
the electron lifetime. 
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I. INTRODUCTION 


The space-dependent spin dynamics in ferromagnetic metals have been the subject of con¬ 
siderable investigation over the past several years. This subject is of great signihcance for 
technological applications such as spin-transfer torque random-access memory (STT-RAM) 
and microwave generators or detectors. In the spintronics devices, which are nanoscale 
magnetic multilayer systems, inhomogeneous magnetization dynamics are caused by an ex¬ 
ternal applied held or thermal excitations due to the hnite size effectsi*^ and the interface 
effects^. The spin torques, which arise in non-equilibrium conditions include some terms 
with different forms^'-, and the space-dependent spin torques caused by inhomogeneous 
magnetization dynamics are different from the current-induced spin torques, namely, the 
spin transfer torque^”— and the spin-orbit torque^^^— . The former torque is due to the 
indirect interaction between magnetization through conduction electrons, while the latter 
torque is due to the direct interaction between magnetization and conduction electrons. 
Generally, the space-dependent torque acting on magnetization at r can be represented by 
J dr'dt M{r,t) x x{r,r',t,t')M{r',t'), where x(r,r',f,f') is the spin susceptibility tensor 
(in detail, see ([9])). We are interested in the dynamic part of this torque, which describes 
the interaction between magnetization M (r) and the time-derivative of other magnetiza¬ 
tion M (r'), and it corresponds to Gilbert damping torque aM x M when magnetization 
processes uniformly, where a is the Gilbert damping constant^. Although a large number 
of studies have considered magnetization damping in uniform precession systems^i^i^— (i.e., 
local damping), little is known about the damping mechanism in inhomogeneous magneti¬ 
zation dynamics2I^i22 (i.e., nonlocal damping). Owing to the development of recent experi¬ 
mental techniques, detailed and unihed theories of nonlocal damping are more desirable for 
spintronics applications. 

Previous theoretical studies show that magnetic damping arising from spin diffusion is 
enhanced by the spin wave that describes the plane wave dynamics of magnetization^^. The 
damping coefficient is represented by a = Oq -f- rjq^, especially in the long wavelength limit, 
where the first term represents damping for uniform precession systems and the second 
term is the contribution from spin wave motion with wave vector q. The coefficient rj is 
the diffusion constant, which depends strongly on the lifetime of the conduction electrons. 
Moreover, experimental evidence for nonlocal damping has been reported by the recent 
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work-, where the Gilbert damping constants for two modes (center mode and edge mode) 


depend on the size of the nanoscale magnets, implying that additional damping torque was 
caused by finite size effects. 

The purpose of our work is to calculate the dynamic part of the spin torque induced by 
the spin wave motion of magnetization in a Rashba spin-orbit coupling (RSOC) system. No 
previous studies have discussed the inhomogeneous magnetization dynamics in SOC systems. 
In our previous work—, we discuss qualitatively the effects of SOC in the RSOC system 
without electrons scattering by impurities. In this paper, we show the results including 
effects of the impurity scattering beyond the previous results in the clean limit and discuss 
the effects of impurity-ladder-sum vertex corrections. 

The reminder of this paper is organized as follows. In Sec. II, we explain our model 
and formulation. In Sec. Ill, we show our calculation results for the dynamic part of the 
spin-wave-induced spin torque. Finally, a summary of our work is given in Sec. IV. 

II. MODEL AND FORMULATION 

The exchange coupling between electrons and magnetization is represented by 



( 1 ) 


where n{r,t) = M {r,t) /M is the unit vector of the time-dependent magnetization at r 


and A is the exchange energy. Here, s (r) = (r) d"ip (r) /2 is the electron spin density 

with the electron held operator (r) = (r), (r)] and Pauli matrix vector <t (the 

hat symbol indicates a 2 x 2 matrix). The direction of magnetization is parallel to 2 ;-axis 
at the equilibrium state, and the contributions from the transverse component (denoted by 


T) are treated perturbatively. Based on linear response theory, the time-dependent energy 
is written as 




where N is the number of unit cells. The spin susceptibility tensor y contains elements 
which are given by 



( 3 ) 
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The effective field is derived from He^{t) = M~^Eex(t) /dn{r,t), and the spin torqne is 
obtained by Tgt {r, t) = —jn {r, t) x iTex (t), where 7 is the gyromagnetic ratio. In the hrst 
approximation, neglecting the prodncts of pertnrbed qnantities, the spin torqne is written 
as 

2A 

Tst (r,t) = 7 — X (s^(r,t)) , (4) 

where = 2 A J dt' J dr'x{r,t,r',t') is the transverse component of the 

electron spin density expectation valne. In slow modnlations of the magnetization motion, 
n-*- (r', t') ^ n-*- (r', t) + (t' — t) (r', t) is satisfied, and the spin torqne is rewritten as 


Tst {r,t) 


dr'n^ 


X 


(3 {r — r') {r', t) + a{r — r') (r', t) 


(5) 


in an electron system with translation symmetry. Here, 

Hr) = x{r,oJ = 0 ), 


( 6 ) 


is the coefficient tensor of the static part, and 


a (r) 


AjA^N 
^ M 


lim , 

(.^->■0 OOJ 


(7) 


is the coefficient tensor of the dynamic part. The Fonrier transformation of spin snsceptibility 
is dehned as x{uj) = J dte'^*x(t). The hrst term on the right-hand side of Eq.(j5]) is the 
RKKY-type interaction, which is regarded as the time-dependent correction for precession 
torqne cansed by the static magnetic held; however, this term is irrelevant for onr work. In 
this stndy, we focns only on the dynamic spin torqne, which is described in second term on 
the right-hand side of Eq.([5]). Expanding the magnetization motion nsing a plane wave with 


wave vector q, the dynamic part of Tgt in qr-space is written as x aqfi^. Here, 

~ 4A2 ^ t ~ t ^ 


( 8 ) 


where hS is the spin angnlar momentnm of magnetization in a nnit volnme v = V/N. 
Using a simple calcnlation, it can be shown that the cj-derivative of the real part of 
Xq = J dry (r) vanishes in a; —)■ 0 . Hereafter, the matrix elements related to the 

2 ;-component of dq are neglected becanse we are not concerned with the longitndinal motion 
in this stndy. 

We calcnlate Xg ( 1 ^) using a Green fnnction techniqne that considers the electron-impurity 
(spin-independent) interaction within the hrst Born approximation^^. The spin susceptibility 
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is obtained using analytic continuation of the Matsubara Green’s function, which is written 


as 


»i/r 


< (■>'„) = / dr,e'-"’> (T4 (r,) s™,) 

Jo 


k n' 


g {iUn', k) (iUn', iWn' + ir'n, q) g (iUJn' + i^n, fc + q) d” 


( 9 ) 


where g{iun,k) = — (T'lpk (g )is the one-particle Green’s function and 

(in;„/, iwn'+ ir'n, q) is the vertex function. This vertex function satisfies the following 
Ward identity: 




P1P2 


(Tipip2 <J2 


iUrf + '\v„, q) (iu„., iWn. + q ), 


( 10 ) 


where 


ff( tiPiP2C’'2 1 i^n' T i^ri) q) 


ep 


TTV-pN 


E 

fc 


5'<7ipi (iWn' + il^n, fc + q) 9p2a2 (W, fc) • (H) 


Here, fc/r = iriy-pnou'^ is the inverse electron lifetime r [up is the density of states per volume 
at the Fermi level, no is the density of impurities, and u is the scattering constant of the 
short range potential uv6{r)]. We define the following 4x4 matrices: 


K = 



^tttt 

nmt 

ntttt 

nmt 

n^tti 

Htm 

Htitt 

Htm 

H^ttt 


n^tt 

ifmt 

V 

n^m 

ifiitt 



Then, the vertex function is given by 


( 12 ) 


0-1(72 


(ia;„/,in;n/ in„, q) = 




UIPIP 20-2 


[ iUn ', iUJn ' + iz^n, Q) C ' 


PlP 2 ■ 


P 1 P 2 


(13) 


Substituting Eq. flT^ into Eq.Q and performing the conventional fc-sum and n'-sum in the 
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low-temperature limit, we obtain the following form of Eq. (l8|) : 


o-q = Q!q(T° -|- + a^(—i)(j^ -I- 


z/pA^ 


2 ^eF 

a = - 

2Se 


= 


F 

z/pA^ 

2S'ep 

, _ Z/fA 2 
2 ^eF 


Re V 

/ ^ \ y y y y / crcrt 

(j 

CT 

im J] (A™ri™ - /?“ri“) _ 

CT 

Re V 

/ ^ \ y y y y / crcri 


(14) 


where (nq^’^'^)aipip2a2 = (cf/ttz/fA) = 0, fc + ^ super¬ 

scripts R and A denote the retarded Green’s function and the advanced Green’s function, 
respectively. 

In Eq. lfT^ . aq is the conventional Gilbert damping coefficient, and a^, which is the 
coefficient of x (—i)(j^n^ = —, is the correction for the equation of motion. Neither 
nor have been examined in previous studies because these terms disappear in non- 
SOG systems and in uniform precession systems. The physical meaning of these terms is 
discussed in the next section. 

The electron system of our model is a 2D electron gas incorporating RSOG^, which is 
described by 

"H = eF ^ V’fc (l^ol^ + Afc ■ (t) V’fc, (15) 

k 

where ep = h‘^kp/2m is the Fermi energy and feg = (kx/kp^ky/kp) is the normalized wave 
vector. Here, is defined as 


Afc = Afco (sin 9k cos 0*,, sin 9k sin 0^, cos 9k ), (16) 

where (j)k = — tan“^ (A;a,/A;j^), 9k = cos“^ (— Ao/Afc(,), and Xk^ = \/ A^fcg -|- Aq. Moreover, 
Aq = and A is the strength of RSOG (the conventional Rashba parameter ap is written 

as Or = Xep/kp). The eigenstate corresponding to the eigenenergy of "H, ek± = cf (^o ^ -^fco)) 
is given by 

|fe, +(-)> = e-<+"^ cos ^ It (t)> + (-)e+l->‘^sin^ |t (t)). (17) 
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Then, the matrix elements of the retarded and advanced Green’s function are written as 


(i^n, k) = Y^ {(^i\a)k k) {aW2)k . 

a=di 

{iujn,k) = {iun +ep - eka±ik/T)~^, (+forR,-for A). (18) 

We can calculate aq from Eq. ffT^ . Eq. lfT^ . Eq. ffTT|) . Eq. (fT7|) . and Eq. lfTSj) . 


III. RESULTS AND DISCUSSION 

In the numerical calculations, we set the following parameters: Aq = 0.1, z/p = l/ep, 
A = 0.3, and S' = 1. We show the numerical results of two cases: h/r = O.Olep and 
Ti/t = O.OSep. However, we do not show the detail results of in this section because we 
have determined that this term is sufficiently small and can be neglected. In our results, the 
maximum absolute value of for Tijr = O.Olep is about 0.07 and that for Ti/t = 0.05ep is 
about 0.06, which satishes <C 1 (i.e., <C n-*-). 


A. Damping torque 


In this subsection, we discuss the results of which is the coefficient of conventional 
damping torque, x n^. The ^-dependence of aq is shown in FigHJ It can be strictly 
shown that aq does not depend on the direction of q. 

In the clean limit. 


Oiq — 


where 


lap (T'o) — 


+ al3\^ {kl + kfj ^ ^ ^ 

&F L (J.2 _ *.2)2 -foF (*) . 

0 [Qo - \K - fc/3|) 0 {K + kp- qo) 
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90 = -r- 
/Cp 


(19) 


( 20 ) 


[(*:„ + k/,) - jJ] [g^ - - kf) ] 

represents the strength of the nesting effects for spin excitations between the a-band and 
the /9-band. Here, ^ 1 + AV2 T + ^V4 is the radius of the Fermi sphere of 

the ±-band. At go = 0, ± A;+, 2/c±, aq diverges because of the strong nesting caused by the 

adjoining Fermi surfaces. However, in a non-SOC system, aq diverges only at go = fc- ± k+— 
because intra-band transitions are forbidden (i.e., the contribution from a = P is zero). The 
divergence at g = 0 in SOC systems is well understood from previous results^^, but the 
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divergence at qo = 2k± arising from the intra-band transitions is a new result related to the 
inhomogeneous dynamics of magnetization. 

From FigiH it is conhrmed that the aq divergence is suppressed by the impurity scattering 
of electrons, and the value at each peak decreases with decreasing electron lifetimes. These 
results imply that the nesting effects are reduced by broadening the Fermi level due to the 
increase in self-energy. There are cases that ag increases with decreasing r because the 
line width of each peak increases with r. These behaviors are also conhrmed in the non- 
SOC system in our previous study2^. This implies that aq increases due to spin diffusion 
originating in inter-band transitions. 

The behaviors of aq near q = 0 are shown in FigI2J The value of at g = 0 corresponding 
to the Gilbert damping constant in a uniform precession system diverges in the clean limit, 
but it converges in the presence of impurities^^. On the other hand, in a non-SOC system, 
aq=o is always zero regardless of the concentrations of impurities; this behavior is required 
from the angular momenta conservation law, which is satished by taking the vertex correction 
(VC). However, in the presence of magnetic (spin-dependent) impurities, 0^=0 has a hnite 
value due to the violation of the conservation law^i^S. in a non-SOC system, the spin- 
dependency of impurities causes dramatic differences in the aq behavior; however, the spin- 
dependence on impurities rarely appears in SOC systems^. Therefore, we do not need to 
provide a detailed analysis of the effects of magnetic impurities in our RSOC system. 

From FigJSl it is clear that results including the VC term are larger than those without 
the VC term, especially at g ~ 0. However, in a non-SOC system, aq with VC is zero at 
g = 0, and aq is smaller than the results without VC. The VC term reduces the overvalue 
of the contributions from spin diffusion originating in inter-band transitions. In an SOC 
system at g ~ 0, the contributions from intra-band transitions are much larger than those 
from inter-band transitions. We conclude that the VC term of an SOC system enhances the 
contributions from nesting effects, which originate in intra-band transitions. 


B. Distorting torque 

In this subsection, we discuss the results of both a^a^ and a^a^. In the RSOC system, 
these terms vanish at g = 0, but not at g > 0. From our calculations, the sum of these terms 





FIG. 1. The g-dependence of Uq (divided by i/pA^/S'er = 10 ^). The results in the clean limit are 
indicated by the dashed line. 



qo 


FIG. 2. Behavior of aq near g = 0. 


is given by 

O^gd “T OigCT — dg ^qi 

where Rq is the reflection matrix in terms of the direction 9q = ia.'n~^{qy/qx)'- 


Rq 


( COS 26q sin 26q 
sin 29q — cos 29q 


( 21 ) 


( 22 ) 
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FIG. 3. Comparison of aq in go < 0-05 including the vertex corrections (VC) and not including 
the vertex corrections (NVC). The solid lines show the results for VC, and the dashed lines show 
the results for NVC. The blue lines are results for Tl/t = O.Olep, and the red lines are results for 
U/t = O.OSep. 


Moreover, is the g-dependent coefficient of new torqne, x Rqh^. The direction of 
this torqne is shown in FigJU We call this the distorting torque because it applies q- 
dependent force to the precession motion. For example, the distorting torque is parallel to 
the damping torque when is perpendicular to g; on the other hand, the distorting torque 
is anti-parallel to the damping torque when is parallel to q. 

The g-dependence of is shown in FigJ3 We confirm that the inequality equation, 
< aq, is satisfied in the entire g-range. Therefore, the precession orbits necessarily 
decrease with time. 

In the clean limit. 


a„ = 


z/pA^ 


A' 


a/3 


{kj-klf hi+ {kl + kfj 
{kl - klf 


lap (go) • 


(23) 


The peaks of a^^ occur at same points with those of aq because the right hand side of Eq. fl2^ 
includes Rp (go)- This equation implies that the distorting torque is unique to the RSOC 
system because a^^ is proportional to the square of A. We find that the contribution from 
the intra-band transitions of a^^, which are given by 2A^r'pA^/c^J„Q(go)/S'ep(fc^ — k\Y, is 
equal to that of aq, but the contribution from the inter-band transitions is not equal. These 
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results reflect the differences in the radius of the Fermi surfaces and in the spin directions 
at the nesting points between the +-band and the —band. These differences are attributed 
to the fc-dependent spin states originating from RSOC. Moreover, the negative values of aq 
at 1.7 ^ go ~ 2 imply that the larger contributions come from inter-band transitions than 
from intra-band transitions. 

In the presence of impurities, the divergence of is suppressed similarly to the results of 
ttq. The behavior of near g = 0 is shown in Figj6l At g = 0, is always zero regardless 
of the presence or absence of impurities. However, in g —)■ 0, diverges (i.e., discontinuity 
at g = 0) in the absence of impurities, but converges to zero (i.e., continuity at g = 0) 
in the presence of impurities. In the absence of impurities, at g < because 

only intra-band transitions occur, but this equality is no longer satished in the absence of 
impurities due to spin diffusion originating in the inter-band transitions. Consequently, 
decreases at the peak near g = 0, and this peak shifts right with decreasing electron lifetimes. 

FigJT] shows the differences between the results including VC and those without VC. The 
values with VC are larger than those without VC, and VC for clearly has effects similar 
to those for aq (see Fig|3]). 



FIG. 4. Direction of x Rqh^. 
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FIG. 5. The g-dependence of (divided by z/pA^/S'eF = 10 ^). The results in the clean limit 
are indicated by the dashed line. 



0.00 0.02 0.04 0.06 0.08 0.10 

qo 

FIG. 6. Behavior of near q = 0. 


C. Oscillation modes 

In the small plane wave field h-^, if we can neglect /3q, then the linearized equation of 
motion is given by 


fig = —yriq X x x . 


(24) 
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VC (0.01 ep) 



0.00 0.01 0.02 0.03 0.04 0.05 




FIG. 7. Comparison of in qq < 0.05 including the vertex corrections (VC) and not including 
the vertex corrections (NVC). The solid lines show the results for VC, and the dashed lines show 
the results for NVC. The blue lines are results for fi/r = O.Olep, and the red lines are results for 
U/t = O.OSep. 


We solve Eq. fl24p by assuming harmonic-time-dependence of h,q(~ without the cor¬ 
rection term . We obtain the eigenequation = Xg , where is the oscillation 

mode corresponding to each component of the magnetic held, (hg is the component of 
right-hand rotation relative to the direction and h~ is the component of the opposing 
direction). Assuming aq,a^^ -C 1, we hnd 

Ug = cos [ojt) ± i sin [ojt) =f sin (ut -|- sin 20^) , (25) 


and 


xt = 


7 


-■ ^ 26 ) 
'yH^ =F ittgCU 

Eq. fl25ll indicates that each magnetization precession describes an elliptical trajectory. 
Stronger SOC, which generates a large magnitude of produces a more elongated elliptical 
orbit. When > 0, the major axis of the elliptical orbit is parallel to q, while the major 
axis is perpendicular to q when < 0. 

From Eq. (|26|l . it is conhrmed that Xq is not affected by when Og, <C 1. Therefore, 
Q!g^ cannot be determined experimentally in the same way as the estimation of ctg, which 
can be measured from the line width of magnetic susceptibility. To experimentally estimate 
the value of direct observation of the elliptical orbit motion is required. 
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We suggest that the time-resolved magneto-optical Ker effect (TRMOKE)^ is an appro¬ 
priate method for detecting the elliptical trajectory in RSOC system. For these observations, 
a large value of is required: both strong RSOC and nesting effects are required. In our 
results for h/r = O.Olep, the maximum value of is about 0.13, which is sufficient for 
detection at g ~ 0.12kp ~ O(10“^)A ^ . However, an artificial excitation of spin wave which 
wave number is larger than 0.01 A ^ is technically difficult at this time. Moreover, for 
clear detection of the qr-dependent precession motion, the spot size of TRMOKE must be 
reduced from 10 pm to a few nanometers, which is currently a difficult target. Thus, new 
experimental techniques are necessary to resolve these issues. 

IV. SUMMARY 

We calculate the dynamic part of the spin torque induced by the plane wave dynamics 
of magnetization in a RSOC system. In addition to the conventional damping torque, 
our results show that distorting torque, which is a phenomenon unique to RSOC systems, 
originates from the inhomogeneous dynamics of magnetization. The magnitudes of these 
torques depend on the strength of the nesting effects, and they are reduced by decreasing 
the electron lifetime. The vertex corrections for these terms are sufficiently large, especially 
in the long wavelength limit, correcting the deficiency of contributions from nesting effects, 
which originate in the intra-band transitions. In the resonant plane wave field, the oscillation 
mode of magnetization precession exhibits an elliptical trajectory whose major axis depends 
on the direction of the wave vector. However, to observe this elliptical precession motion, 
improvements in measurement sensitivity and experimental techniques are required. 
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